%% Comparison of Errors: 

JLF = LJMM_F;
JLS = LJMM_SC;
errJLF = reshape(JLF./u_vec - 1,[N,T]);
errJLS = reshape(JLS./u_vec - 1,[N,T]);

LJF = max(abs((errJLF(:,T))));
LJS = max(abs((errJLS(:,T))));

LJF2 = norm((errJLF(:,T)),2)/sqrt(N);
LJS2 = norm((errJLS(:,T)),2)/sqrt(N);

ErrorJL = [LJF, LJS,LJF2, LJS2];
